Power Spectrum Precision for Redshift Space Distortions 
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Redshift space distortions in galaxy clustering offer a promising technique for probing the growth 
rate of structure and testing dark energy properties and gravity. We consider the issue of to what 
accuracy they need to be modeled in order not to unduly bias cosmological conclusions. Fitting for 
nonlinear and redshift space corrections to the linear theory real space density power spectrum in 
bins in wavemode, we analyze both the effect of marginalizing over these corrections and of the bias 
due to not correcting them fully. While naively subpercent accuracy is required to avoid bias in 
the fixed case, in the fitting approach the Kwan-Lewis-Linder reconstruction function for redshift 
space distortions is found to be accurately selfcalibrated with little degradation in dark energy and 
gravity parameter estimation for a next generation galaxy redshift survey such as BigBOSS. 



I. INTRODUCTION 

The pattern of galaxy clustering in three dimensions, 
and its evolution, encodes abundant information on the 
cosmological parameters affecting matter growth. Ongo- 
ing and next generation spectroscopic galaxy surveys will 
vastly increase our measurements of this clustering, and 
our knowledge of cosmology if we can accurately inter- 
pret the results in terms of theory. Measurements accrue 
an extra contribution to the redshift, and hence appar- 
ent position along the sight, from the galaxy peculiar 
velocities induced by the inhomogeneous density field; 
this gives rise to an anisotropy in the observed clustering 
known as redshift space distortions (RSD). 

These distortions carry information on the growth rate, 
as opposed to just the growth amplitude, and so are valu- 
able for probing cosmology, as well as the gravitational 
strength driving the growth. However, linear theory is 
insufficient for accurate relation of the redshift space 
galaxy power spectrum to the true (real space) matter 
density power spectrum, even on quite large scales, or 
wavenumbers k > 0.05 ft./Mpc, where the vast major- 
ity of the statistical leverage lies [1-6]. Numerous cor- 
rections involving higher order perturbation theory have 
been employed [7-10] that extend the validity but the 
region k > 0.1/i/Mpc is still problematic, especially for 
quantities involving the growth rate and the gravitational 
growth characterization. For example, [6] demonstrates 
that these first principles approaches deliver results for 
the growth rate that are biased by several standard de- 
viations when using modes out to k = 0.2 h/Mpc. 

Here we investigate a basic question: how accurately 
does one actually need to know the redshift space dis- 
tortion mapping in order to extract the cosmological and 
gravitational parameter information without substantial 
bias or degradation? Similar questions have been consid- 
ered for weak gravitational lensing, for example, where 
one asks how well the nonlinear matter power spectrum 
needs to be known to estimate cosmology from the lens- 



ing shear power spectrum [11, 12]. 

In Section II we introduce the correction, or recon- 
struction, function for the redshift space power spectrum 
and review the KLL [6] form for it. Section III uses the 
Fisher bias method to compute both the individual pa- 
rameter bias and joint confidence contour bias due to 
misestimated RSD, thus giving criteria for the accuracy 
to which the RSD effects must be known. Adding fit pa- 
rameters for uncertainties in the reconstruction function 
in Sec. IV, we assess the impact of marginalizing over 
them on the cosmological parameters, in particular for 
tests of dark energy and gravity. Section V summarizes 
the results and conclusions. 



II. GALAXY POWER SPECTRUM 

A. Redshift Space Power Spectrum 

In real space the matter density power spectrum is ex- 
pected to be isotropic, and the linear power spectrum 
grows in a scale independent manner through the growth 
factor D(z), where z is the redshift. The observed, red- 
shift space galaxy power spectrum involves a transforma- 
tion to redshift space due to the velocity effects, and a 
bias relation b(z), usually taken to be scale independent, 
converting the dark matter overdensity to galaxy over- 
density, and the effects of nonlinear structure formation. 
Each of these is modeled in various ways, with attendant 
uncertainties. 

We write the anisotropic redshift space galaxy power 
spectrum as 

P{k, fj,, z) = P r {k, z) M(k, fi, z) F(k, fi, z) (1) 

where P r is the isotropic real space matter power spec- 
trum, M is an approximate model for redshift space dis- 
tortions (including galaxy bias), and F is the reconstruc- 
tion function accounting for nonlincarities and more ex- 
act velocity effects. 
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The linear mass power spectrum P r is given by a Boltz- 
mann numerical code such as CAMB [13]. It depends 
on the cosmological parameters through its shape as a 
function of wavenumber k and through the growth factor 
D(z) giving its amplitude evolution. Since we will cor- 
rect the RSD modeling by the reconstruction function, 
we choose M to be simply given by the linear theory 
prediction, the Kaiser approximation [14], 



M{k,n,z)= [b(z) + f(z)n 



2l2 



(2) 



where b is the galaxy bias, / = d\nD/d\na the growth 
rate of density perturbations that in the linear regime 
grow as 5 ~ D(a), where the scale factor a = 1/(1 + z), 
and n is the cosine of the angle made by the perturbation 
wavevector k with respect to the line of sight. Beyond 
the linear regime, b could be scale dependent, i.e. b(k), 
but we will absorb that possibility into the reconstruction 
function. The reconstruction function is fitted to N-body 
simulations by the analytic form of Kwan, Lewis, & Lin- 
der (KLL; [6]), 



F(k,^z) = 



A(k,z) 



1 + B{k,z)k 2 fj, 2 



C(k,z)k 2 fi 2 



(3) 



This form has been found to reproduce accurately re- 
sults of N-body simulations over a wide range of redshifts, 
and for halos of various masses as well as dark matter; see 
[6, 15] for details. Note that A, B, C may be cosmology 
dependent, just as / and P r are, and their universal- 
ity is a subject of ongoing research, but here we treat 
them as independent parameters (as an analogy, recall 
how people treat coefficients within Halofit also as uni- 
versal, though here we let the values of A, B, C float; 
also see Sec. IV D). The factor A characterizes nonlinear- 
ity of the real space power spectrum, B describes velocity 
effects such as damping from a Lorcntzian velocity dis- 
persion but also higher order multipole terms, while C 
describes nonlinear enhancement for large k/j, and possi- 
bly breaks the degeneracy in the two roles of B. 



B. Galaxy Clustering Information 

The cosmological information inherent in the galaxy 
power spectrum can be estimated through the Fisher in- 
formation matrix. The full set of parameters {pi} in- 
cludes the cosmological parameters, astrophysical param- 
eters such as galaxy bias, and parameters for the re- 
construction function. Sensitivity to cosmology enters 
through the derivatives dP/dpi and the error covariance 
matrix for the redshift space galaxy power spectrum P. 

We follow the usual prescription [16-18] where the co- 
variance matrix comes from sample variance (finite vol- 
ume) and shot noise (finite resolution of the density field 
by sparse galaxies). Taken together, the error can be 
thought of as depending on the number of modes that 
the galaxy redshift survey samples. Treated as Poisson 



sampling of the density field, the statistical error is 

( jp = P + n- 1 (4) 

from these two effects. The number of Fourier modes is 
the volume of a Fourier cell times the number of cells, 

^modes = 2nk 2 dk dfi x V s 

urvcy 

/(27T) 3 . (5) 

Therefore the error covariance matrix C is 

(6) 



C = P 2 



1 + nP 



nP 



8tt 2 



Vsurveyk dkdf-t 



''survey r 

Since the Fisher information matrix is constructed 
from C _1 multiplied by the sensitivity derivatives 
dP/dpi, we can use the P 2 factor to convert the deriva- 
tives to involve In P, which will be useful in treating the 
multiplicative factors in Eq. (1). In summary, the Fisher 
matrix is 

k 2 AkA/j, 



8tt 2 



(7) 

where the survey volume is reduced by the shot noise to 
a z, k, and (i dependent effective volume 

2 



VeB(k,fJl,z) = V surV cy(z) 



n(z)P(k, fi, z) 



n(z)P(k, fj,, z) 



1 



(8) 



When the galaxies densely sample the underlying field, 
the effective volume approaches the survey volume, oth- 
erwise modes are lost, diluting the effective volume due 
to increased noise. 

Note that the logarithmic derivatives can be written 

as 

dlnP dlnP _ (<91nP r + dlnM + dbxF) 



dpi dpj dpi 

{d\nP r + dhxM + dhxF) 



1 



so only the <91nP term depends on A, B, and C. 
The reconstruction function derivatives are 
dF 1 
dA 
dF 
dB 
dF 
dC 



Bk 2 ^ 2 
-Ak 2 fi 2 



{1 + Bk 2 n 2 ) 2 
k 2 n 2 , 



(9) 

(10) 

(11) 
(12) 



and are otherwise taken not to depend on cosmology. 
This is because we use A, B, C purely as fiducial values, 
and investigate how their variation (from astrophysics or 
cosmology) impacts the cosmological parameter estima- 
tion. 

Our fiducial case attempts to match F to the simula- 
tion results in [6], with estimated 

4 1.58 

A(k) = 1 



B(k) 
C(k) 



0.39 h/Mpc, 
20 (Mpc/h) 2 

8e -fc/(0.176VMpc) ( Mpc//l )2 



(13) 

(14) 
(15) 
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The resulting redshift space distortion reconstruction 
function of Eq. (3) is shown in Fig. 1. We emphasize 
that these are merely the fiducials; we allow the values of 
A, B, C to float freely in bins of wavenumber. This pro- 
vides a model independent variation of the power spec- 
trum (within the reconstruction form) and we can then 
investigate the influence of such variations on the cos- 
mological parameter estimation. Conversely, the ques- 
tion can be phrased as "what is the accuracy required 
on knowledge of the galaxy power spectrum in order to 
deduce the cosmology with confidence?" 

We later contrast this fiducial with fiducial (A, B, C) = 
(1,0,0), i.e. assuming that perturbation theory (for ex- 
ample linear theory in the Kaiser case, although F also 
works with higher order perturbation theory [6]) fully 
captures RSD effects in the model M. 
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mately 3 x 10 -4 h 3 Mpc~ 3 . For the exact distribution 
adopted see Table I. There are actually two populations 
of galaxies: luminous red galaxies (LRG) and emission 
line galaxies (ELG), each with their own galaxy bias 
value. These galaxy biases are taken as free parameters 
to be marginalized over, for each redshift bin of width 
0.1. Their fiducials are b(z) = b D(z = 0)/D(z) with 
frELG _ o.8 and b}; RG = 1.6, which provide good fits to 
observations. Galaxy populations with different biases 
can help reduce sample variance [20], with the Fisher 
matrix involving a sum over populations, i.e. 



E 

XY 



d\nP x dlnPy 

dpi dpj 



n x Px 



n x P. 



x 



1 



nyPy 



UyPy + 1 



(16) 



Note that for multiple populations the factor V c r in 
Eq. (8) involves the shot noise, i.e. nP 1 of each popu- 
lation. 
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22.6 
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8.45 


3.04 


0.35 


4.02 


3.07 


0.45 


2.65 


3.09 


0.55 


2.99 


3.10 


0.65 


3.99 


3.11 


0.75 


5.15 


3.12 


0.85 


5.36 


1.89 


0.95 


5.02 


0.33 


1.05 


4.80 


0.04 


1.15 


4.49 




1.25 


4.04 




1.35 


3.02 




1.45 


2.00 




1.55 


1.15 




1.65 


0.43 




1.75 


0.12 





TABLE I. Spectroscopic survey number densities adopted for 
emission line galaxies and luminous red galaxies, in units of 
10" 4 ft 3 /Mpc 3 , for each redshift shell. 



FIG. 1. The redshift space distortion reconstruction function 
F(k,n) is plotted for the fiducial expressions for A, B, C for 
three values of angle /j,. 

The analysis is carried out in the next sections in two 
ways: in Sec. Ill we compute the effect that a given level 
of unrecognized power spectrum deviation in some k bin, 
i.e. a systematic error in modeling, has in biasing the 
cosmological conclusions, and in Sec. IV we recognize the 
existence of systematic uncertainties and treat them by 
marginalizing over the A, B, C values for each k bin. 

C. Survey Characteristics and Parameters 

For the galaxy redshift survey data we consider a 
next generation spectroscopic survey of the quality pro- 
posed for BigBOSS [19], covering 14000 deg 2 from z — 
0.1 — 1.8, with a galaxy number density n of approxi- 



For the cosmological parameters we use the physical 
baryon density Q, h 2 and physical cold dark matter den- 
sity tt c h 2 , reduced Hubble constant h, scalar perturba- 
tion tilt n s and amplitude A Sl dark energy equation of 
state parameters wq and w a , and gravitational growth in- 
dex 7. The gravitational growth index gives an accurate 
description of the growth rate for both general relativity 
and a range of modified gravity models, and looking for 
deviations from its general relativistic value of 7 = 0.55 
acts as a test of gravity [21, 22]. The growth index is 
treated as an independent parameter (not a function of 
u>o, w a ) and determines the growth factor at scale factor 
a = l/(l + z), 

D(a) = e fo( da '/ a ')n m (a'r ^ ( 17 ) 

that in this ansatz is used to convert the linear power 
spectrum delivered by CAMB at z = to another red- 
shift z, to account for the effects of modified gravitational 
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growth. Note that the growth rate / = r2 m (a) 7 , and red- 
shift space distortions were highlighted as a test of gravity 
in [23]. 

For the central question of RSD uncertainties we em- 
ploy up to 12 free parameters, taking A, B, C with inde- 
pendent values in each bin of width 0.1 in wavenumber 
above k = 0.1 h/Mpc out to some /c max . This corresponds 
to uncertainties APfc. For the current work we follow 
[11, 12] and consider the uncertainties only as a func- 
tion of wavenumber, not redshift, except in Sec. IV C; 
we also take the KLL form to be accurate while allowing 
freedom in the parameters A, B, C . In summary we fit 
for 8 cosmological parameters and up to 39 systematics 
parameters. 



III. FISHER BIAS 

The first question we are interested in answering is 
what is the sensitivity of the cosmological parameter de- 
termination to errors in modeling RSD. One might have 
M or F wrong, but if this does not mimic a change in 
cosmology then no harm is done. The Fisher bias for- 
malism (see, e.g., [24, 25]) propagates misestimation of 
the observable or theoretical prediction, in this case the 
redshift space power spectrum, into biases on the fit pa- 
rameters. Specifically, we consider the effect of errors in 
the k bin values of A, B, C. 

The Fisher bias on a parameter pi from misestimating 
parameter p a is 



5 Pl = S Pa > (F ) ij (F tu11 ) 



j 



(18) 



where Sp a = p a (truc) — p a (fiducial). The superscript 
"sub" denotes the Fisher submatrix without entries for 
the parameters whose misestimation we are studying. 
(For the specific case here, the submatrix will be 35 x 35 
for the cosmology and galaxy bias parameters, and the 
full matrix adds the reconstruction parameters one at a 
time. We later consider all the reconstruction parameters 
at once.) By evaluating the ratio dpi/dp a for a = A, B,C 
we can assess the sensitivity of the parameter estimation 
to the RSD modeling. To take a weak lensing example, 
[26] found that a particular form of matter power spec- 
trum distortion with amplitude Anl at high k distorted 
estimation of w a derived from shear power spectrum mea- 
surement by a leverage factor of 18: a misestimation of 
10% in Anl yielded a 1.8a bias in w a . 

The bias S P i can be compared to the statistical uncer- 
tainty cr(pi) on the parameter, either directly or through 
the risk statistic 



Ri = Ja 2 (p i ) + 5p 2 i 



(19) 



Treating the bias as a systematic error in this way, to 
restrict the degradation in the statistical error to less 
than 20%, say, requires Spi/a(pi) < 0.66, thus putting a 




constraint on the allowable modeling error 5p a etc. We 
examine two, converse statistics: the cosmological degra- 
dation caused by a certain fractional misestimation of the 
reconstruction parameters Sp a /p a , and the requirement 
on the reconstruction parameter to bound the cosmo- 
logical parameter bias to less than a given factor of the 
statistical uncertainty, Spijoi. These are respectively 



(20) 
(21) 



Figure 2, left panel, shows the matrix of degradations 
Ri/(Ji for fixed Sp a /p a = 0.01 (i.e. 1% uncertainty on 
the reconstruction parameters), where the columns are 
the dark cosmological parameters and the rows are the 
reconstruction parameters. The right panel gives a simi- 
lar matrix of the reconstruction requirements Sp a /p a for 
fixed 5pi/(Ji = 1 (which corresponds to Rijoi = 1.41). 
One can scale the results for different fixed values accord- 
ing to the above equations. The stripe structure arises 
because the bias from A .45, where the subscript indi- 
cates the center of the k bin, is of opposite sign from 
that of j4o.25j and A0.35 lies in between near null effect, 
and similar for B and C. 

The degradations in determination of the dark param- 
eters wo, w a , 7 are less than 22% for 1% shifts in the re- 
construction parameters in all cases except A). 25, A).45j 
and -Bo. 45- For A0.25 and A). 45 the risk error on w a can 
exceed the statistical uncertainty by a factor 3. For the B 
parameters the worst case is degradation by 1.5. These 
results offer indications of what physics must be most 
accurately understood, i.e. the nonlinearity from A and, 
somewhat less critically, the velocity effects from B. 

In the converse analysis of what accuracy is required 
on the reconstruction parameters to ensure that a bias 
is restricted to below la, we find that 5% accuracy is 
sufficient for all parameters except for the A^, plus Bo. 25 
and i?o. 45- Knowledge of A0.25 and A0.45 are needed to 
0.3%, £0.45 to 0.9%, £0.25 to 1.4%, A).35 to 1.5%, and 
A0.15 to 3.1%. 

While these approaches give indications of sensitivity, 
they treat the cosmological parameters one by one while 
a power spectrum misestimation will generally impact 
several at once. This can either tighten or loosen over- 
all requirements, depending on the covariances. To take 
this into account we use the A% 2 method [27, 28]. This 
describes the fuller impact of biasing cosmology through 
quantifying how far from the fiducial the best fit cosmol- 
ogy is shifted relative to the confidence contour, taking 
into account degeneracies between the reconstruction and 
cosmological parameters. This measure is given by 



(22) 



where the sum runs only over the reduced parameter set 
whose bias we are interested in, e.g. wq and w a for a 2D 
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FIG. 2. [Left panel] The ratio of the root mean squared error, 
or risk, to the statistical uncertainty, Ri/ai, is plotted for each 
dark cosmological parameter in the case of a 1% deviation 
in a reconstruction parameter. [Right panel] The fractional 
requirement on each reconstruction parameter 5p a /p a needed 
to ensure bias less than la, i.e. 8pi/oi < 1 is plotted. Dark red 
indicates danger (high risk or tight requirement), with lighter 
colors showing reduced impact. For the left panel the color 
scale is Ri/cn > 2 (dark red), 1.4-2 (medium orange), 1.05- 
1.4 (light yellow), and 1-1.05 (white). The right panel has 
\Sp a /pa\ < 0.01 (dark red), 0.01-0.05 (medium orange), 0.05- 
0.2 (light yellow), > 0.2 (white). Here fc max = 0.5/i/Mpc. 



wo-w a joint likelihood contour plot. The reduced Fisher 
matrix 

F (red) 

is marginalized over all other cosmological 
and galaxy bias parameters (the reconstruction param- 
eters have already been taken into account in obtaining 
S wq etc.). The bias A% 2 accounts for the property that 
biases in, say, the direction of the narrow axis of the 
Fisher ellipse are more detrimental than those along the 
degeneracy direction. 

Figure 3 illustrates the 2D bias induced in the dark 
energy and growth parameters, here for a 1% misestima- 
tion in the reconstruction parameters one by one. Most 
such reconstruction parameter errors do not significantly 
affect the joint parameter likelihood. In the WQ-w a plane, 
none of the C parameters and three of the B parameters 
do not bias the best fit outside the la contour, and -Bo. 45 
remains within the 2a contour. Only errors on the A0.25 
and A0.45 parameters are particularly damaging, causing 
a bias of up to Ax 2 = 39 (approximately at the 6a level). 
The covariancc between the shifts in wq and w a is crucial; 
if the same bias in w a and an even larger bias in wo oc- 
curred such that the joint shift lay along the degeneracy 



axis, then the 2D bias would be scarcely outside the 2a 
uncertainty contour. For the w a -"f plane, the biases are 
less severe, with only A0.25 and A0.45 causing more than 
a 2a joint bias, at A% 2 ~ 16. 

Treating the reconstruction errors one by one effec- 
tively takes a localized bump in the reconstruction func- 
tion. A smooth variation would instead affect several 
of the reconstruction parameters at once; this has a dif- 
ferent effect on the cosmological parameter bias. As an 
example, we simultaneously vary all four A parameters 
by 1%. Since j4o.25 and A0.45 have nearly opposite ef- 
fects this reduces substantially the A% 2 due to varying 
just one of them, e.g. from 39 to 7.7. The 2D bias due 
to such smooth variation is shown in the figures by the 
magenta squares. 

While we have thus far been model independent in tak- 
ing A, B, C to be independent from one k bin to the next, 
we can now consider the difference between two fiducial 
models for the overall reconstruction function F. This 
then includes the effects of variations at all k's simul- 
taneously, and allows a study of bias as a function of 
fcmax- As we consider each successive bin at higher fc, 
we increase the number of modes, reducing the statisti- 
cal uncertainty, but also often increase the deviation in 
the power spectra, increasing the bias in the cosmological 
parameters if we assume the wrong fiducial as the truth. 
The truth is taken to be F as given by Eqs. (13)-(15) 
in Eq. (3), while the incorrect assumption is pure lin- 
ear theory, i.e. simply the Kaiser form for rcdshift space 
distortions, equivalent to A = 1, B = C = 0. 

This miscstimation of the rcdshift space galaxy power 
spectrum causes a bias in cosmological parameters of 



5pi 



F 



sub\ 



' 1 J 

dlnP 

dp. 



P(A,B,C) -P(1,0,0) 



z k 



P(A,B,C) 



k 2 AkAfi 
~~8^ 



(23) 



The systematic biases tend to worsen with increasing 
fcmax, reaching 1.4 in wq, —8 in w a , and 0.2 in 7 for 
fcmax = 0.5 /i/Mpc, and are much larger than the statisti- 
cal uncertainties for all fc max . This is not surprising since 
-Fkaiser can deviate by a factor 2 from the KLL form. 
Thus, neglecting the uncertainties in the reconstruction 
parameters is not a viable option: we must take them 
into account. 



IV. MARGINALIZATION AND 
SELFCALIBRATION 

As an alternative to requiring the power spectrum to 
subpercent accuracy and computing the bias from mises- 
timated reconstruction parameters, we can fit for those 
parameters and calculate the increased uncertainty in 
cosmological parameters due to marginalization over the 
extra inputs. The basic question is how well the model 
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FIG. 3. The biases in the wo~w a and w a -~f planes due to 1% 
misestimation in the 12 reconstruction parameters, one by 
one, are shown by x's (along with the A% 2 if larger than 2.8). 
The contours give the joint 2D la and 2a confidence levels 
on the cosmological parameters when the reconstruction pa- 
rameters take their fiducial ( "true" ) values. Magenta squares 
show the biases when varying all bins of A simultaneously; 
such smooth variations are much less damaging, e.g. reducing 
the individual A\ 2 = 39 and 34 biases to a joint Ax 2 = 7.7 
offset (or the 16 and 15 in the w a -y panel to 2.8). 



needs to be known for precision determination of cos- 
mology with RSD. This is similar to what [11, 12] did 
for matter power spectrum uncertainties applied to weak 
lcnsing cosmology. They used fractional power spectrum 
uncertainties in wavemmiber bins, assumed constant with 
redshift, and applied some level of priors. 



A. Global Fit 

Now our quantities A/-, Bk, Ck in each wavenumbcr 
bin become fit parameters. Again, we can study the ef- 
fects as we extend fc max , using more bins and hence more 
parameters. Including these parameters means that we 
will not be biased any more with respect to the fiducial, 
but the enlarged parameter space will lead to some level 
of degradation of the statistical uncertainties, relative to 
fixing the reconstruction parameters, at the same fc max . 

Table II shows the effect of extending the data to 
higher fc max , while simultaneously allowing for the ad- 
ditional reconstruction parameters in each bin. Despite 
the additional degrees of freedom in the fit, the cosmolog- 
ical parameter estimation sharpens with increasing A; max . 
As long as the form of the reconstruction function holds, 
we obtain an accurate and unbiased cosmology even al- 
lowing for fitting variation in the amplitudes of A, B, 
C in each k bin. This is an extremely promising initial 
result for use of the reconstruction. 

To better understand why the added fit parameters do 
not cause an overall degradation, we look at the corre- 
lation matrix of the 47 parameters in Fig. 4. The block 
of parameters 36-47, representing the reconstruction pa- 
rameters, is not highly correlated with other parameters: 
correlation coefficients are below 0.58 (0.38 for parame- 
ters other than n s ). (Even within the block, only -Bo. 15 
and Co. 15, other than between the Ak, have a correlation 
coefficient exceeding 0.8, reaching 0.90; this is expected 
since for a low k expansion both B and C contribute 
as [i 2 .) This means that the change in power spectrum 
shape due to adjusting the amplitudes of these parame- 
ters in F is not degenerate with a change due to wq or 
other such parameters. That is, the influence of these pa- 
rameters have different k and fi dependences than those 
of cosmological parameters and so we find they can be 
separately fit. 

Moreover, the reconstruction parameters are selfcali- 
brated by the data to good precision. All are determined 
to better than 3% (except Co. 15, to 8%) and most to sub- 
percent level. These propagate into determination of the 
power spectrum to the subpercent level for variation of 
each one individually by ler, except for the extreme cases 
of fj, = 1 and Bo. 15 (Co. 15) which gives 1.1% (1.2%) un- 
certainty. Most combinations, however, give subpercent 
precision. Adding all their uncertainties in the most un- 
favorable way generates an extreme of 2.6% power spec- 
trum uncertainty. Thus unlike the weak lensing probe 
analyzed by [11, 12], redshift space distortions do not 
require any priors to be placed on the power spectrum 
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0.00141 


0.00478 


0.0142 


0.0119 
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0.306 
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TABLE II. la constraints from future galaxy power spectrum data on cosmological parameters, marginalized over galaxy bias 
and redshift space distortion reconstruction. Note Q, m is a derived parameter; fc max is in units of /i/Mpc. Despite the addition 
of more fit parameters when increasing fc max , the cosmological parameters can be better determined. 
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FIG. 4. Correlation matrix of the 47 parameters for fc max = 
0.5 /i/Mpc is shown with color shading reflecting the absolute 
value of the correlation coefficient rij = Cy / \l CuCjj . The 
correlation matrix is mostly block diagonal and the cosmo- 
logical parameters are not strongly correlated with the recon- 
struction (or galaxy bias) parameters, so marginalization does 
not badly degrade cosmological parameter estimation. 



parameters (assuming that the KLL reconstruction form 
is valid). 

Remarkably, in addition to selfcalibration, the addi- 
tional fit parameters have little impact on the cosmolog- 
ical parameter estimation, enlarging the uncertainties by 
only 9%, 22%, and 7% on wq, w a , and 7. And of course 
including the extra parameters removes any cosmology 
bias as suffered in the previous section (modulo model 
validity). Regarding the 12 extra reconstruction param- 
eters, from Fig. 1 we see that F is smooth in k so taking 
bins of width 0.1 in k is reasonable. For completeness, 
for bins of width 0.02 (and hence 60 extra parameters) 
we find that cosmological parameter estimation is mildly 
degraded, with uncertainties ontoo, w a , 7 increasing rel- 
ative to 0.1 width by 16%, 33%, 17%. 



Maximum Wavenumber 



Let us examine the dependence of the results on the 



maximum wavenumber fc n 



used. Note that for the 



fcmax = 0.1 case, the cosmology parameters are not de- 
termined particularly well even though no reconstruction 
parameters are used for k < 0.1/i/Mpc. This is due to 
strong covariance with the 27 galaxy bias parameters. 
Once beyond the linear regime, this degeneracy is broken 
and the correlation coefficients drop, greatly improving 
the cosmological parameter determination (e.g. by fac- 
tors of 2.9-3.4 on the dark parameters, for fc max = 0.2 
relative to fc max = 0.1). This continues for higher fc max , 
despite the addition of further reconstruction parameters, 
but gradually saturates. For example, while relative to 
the fc max = 0.5 case the uncertainties on wq, w a , or 7 at 
tax = 0.2 are greater by a factor ~ 3, at fcmax = 0.3 
the factor is 1.6, and at fc max = 0.4 is 1.15, as seen in 
Fig. 5. Thus, having an accurate reconstruction form 
out to fcmax ~ 0.4 — 0.5 is sufficient for robust cosmolog- 
ical parameter estimation, while selfcalibration obviates 
the need for any prior knowledge of the values of the 
reconstruction parameters. 

Binning such as we use is model independent and 
closest to the weak lensing work. This model indepen- 
dence is important since in the absence of a large suite 
of simulations we may have no particular confidence in 
paramctrizations such as those in Eqs. (13-15). Recall 
that those equations merely give the fiducial values in 
each k bin, and then we allow the bin values to float 
freely and marginalize over them. Given simulations we 
might be able to adopt specific forms and fit for a reduced 
set of parameters, e.g. the coefficients in those equations. 



C. Redshift Dependence 

To give a first indication of whether adding redshift de- 
pendence to A, B, C affects the conclusions we include 
a variation of the characteristic wavenumber scale enter- 
ing in the nonlinearity amplitude A in Eq. (13), i.e. the 
0.39 h/Mpc, writing this as 



fc*(z) = 0.39(1 + z) 



a/1.58 



ft/Mpc , 



(24) 



and adding this evolution parameter cv to the fit. The 
simulation results in [15] indicate that relative to A, the 
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FIG. 5. Extending fc max to values above 0.1 /i/Mpc breaks de- 
generacies, leading to improvements in cosmological parame- 
ter estimation as shown here, even given the addition of recon- 
struction parameters to marginalize over. Reconstruction to 
fcmax = 0.4 — 0.5 ft/Mpc is sufficient to plateau the cosmology 
estimation precision. 



parameters B and C have negligible additional rcdshift 
dependence. Therefore we scale B and C by the same 
factor as A, i.e. 



0.58 



B{k, z) = B{k, 0) A(k, z) /A(k, 0) 



(25) 



B(t|0) 1 + W M)- S ( 1 + ,)- i(M1 



and the same for C. 

The introduction of redshift dependence through a has 
little impact on the cosmological parameter estimation; 
the largest correlation coefficient of a cosmologically is 
0.31, with 7, and overall 0.83 with I?o.45> while a itself is 
determined to within 0.025. Figure 6 shows the influence 
on dark cosmology parameter estimation of marginaliza- 
tion over the reconstruction parameters with and with- 
out rcdshift dependence, and fixing the reconstruction 
parameters (i.e. with a total of 48, 47, or 35 parameters). 
Uncertainties on wo, w a , and 7 increase by only 0.8%, 
0.2%, 5% respectively upon including a. Other forms 
of redshift dependence may have different cosmological 
impact, and this deserves further analysis through simu- 
lations, but the scaling of the characteristic wavenumber 
as used here should give a reasonable first indication. 



0.56 
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0.52 



1 1 1 1 1 1 1 1 r 

- - all marg 

— fixed a, marg F 
fixed F 




-0.4 



-0.2 



0.2 



0.4 



FIG. 6. Joint 2D la confidence contours on the dark cosmol- 
ogy parameters are shown for the cases of all reconstruction 
parameters being fixed (dotted red), marginalized over with- 
out redshift dependence as in most of the article (solid black), 
and additionally marginalizing over a redshift evolution pa- 
rameter a (dashed blue). Here fc max = 0.5/i/Mpc. Note that 
the fixed F case, shown here centered on the true cosmology, 
could be strongly biased if F was misestimated (see Sec. III). 



D. Nonlinear Power Spectrum 

The greatest effect of uncertainty in the reconstruction 
function comes from the parameters A^, as seen in Fig. 1 
and in Sec. Ill regarding the parameter bias. Recall that 
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A(k) arises from the nonlinear effects in the density field, 
and even exists when /i = 0. In this limit A{k) acts to 
map the linear real space density power spectrum to the 
nonlinear regime. Therefore, if we had a robust nonlinear 
(or quasilinear) real space power spectrum then we would 
have no need of a separate parameter as then A(k) = I 
(this has been tested and found accurate to subpercent 
level by [15]). Substantial effort is going into developing 
cosmic emulators [29] that could provide accurate non- 
linear power spectra, eventually including the full set of 
cosmological parameters and redshifts considered here. 
Since that is still in the future, we consider the nonlin- 
ear prescription of Halofit [30] to get an indication of the 
potential impact on our conclusions. 

The linear power spectrum at a given redshift is fed 
into Halofit to give the approximate nonlinear form. This 
removes the need for A(k), setting this equal to one for 
all k and z. As indicated earlier in this section, the simu- 
lation results from [15] imply that for such a normalized 
A then the quantities B and C become substantially red- 
shift independent. Therefore we do not need any hypo- 
thetical model such as the a parametrization, making 
the entire analysis more robust. Furthermore, Halofit 
includes cosmology dependence and so no assumption 
about universality of A is needed. 

We show the results for cosmological parameter esti- 
mation in Table III, for k max = 0.5 h/Mpc, for the three 
cases of using the model independent approach of fitting 
for A(k) in bins, using Halofit and A = 1, and using the 
revised version of Halofit from [31] and A = 1. In all 
cases we still fit for the binned values of B and C. The 
use of functional forms for the nonlinear real space power 
spectrum allows better determination of the cosmological 
parameters, by 12%, 28%, 12% for wq, w a , 7 respectively 
(15%, 33%, 15% for revised Halofit, which has slightly 
more quasilinear power). This offers some promise for 
the future use of cosmic emulators, but in this paper we 
prefer to be conservative in the estimations and use the 
model independent approach. 





Wo 


W a 


7 


\l-rn 


Fit A(k) 


0.0713 


0.306 


0.0163 


0.00794 


Halofit 


0.0624 


0.220 


0.0143 


0.00678 


New Halofit 


0.0603 


0.206 


0.0139 


0.00608 



TABLE III. la constraints as in Table II, using fc ma x = 
0.5/i/Mpc, but for three different methods of treating non- 
linearities. 



V. CONCLUSIONS 

With the ability to map galaxy clustering in three di- 
mensions over large volumes of the universe comes the 
necessity for accurate theoretical interpretation. This en- 
tails linking the isotropic, linear theory real space den- 
sity power spectrum to the observed anisotropic, non- 
linear redshift space galaxy power spectrum. We have 



investigated here some of the relevant issues involving 
nonlinearities in the density field and velocity effects, us- 
ing the Kwan-Lewis-Lindcr analytic redshift space recon- 
struction function calibrated from numerical simulations. 

The main question addressed is to what accuracy 
the anisotropic redshift space power spectrum must be 
known in order to achieve robust cosmological conclu- 
sions. We propagate uncertainties in the power spectrum 
through a model independent binning of reconstruction 
amplitudes with wavenumber and assess the effects of 
deviations from fiducial values. To avoid biasing cos- 
mological parameters such as the dark energy equation 
of state and gravitational growth index requires down 
to 0.3% accuracy on the reconstruction parameters in 
the most stringent cases, while smoother deviations give 
more tractable requirements. Note that it is only those 
deviations that mimic cosmological variations that are 
most important. 

A more flexible and robust approach is to carry out 
a global fit for the binned reconstruction parameters si- 
multaneously with the cosmological parameters, which 
avoids biasing the results so long as the form of the KLL 
function is a good approximation. With 8 cosmology pa- 
rameters and 39 systematics parameters we find that a 
next generation galaxy redshift survey such as BigBOSS 
can tightly and accurately constrain cosmology, for ex- 
ample determining the equation of state time variation 
w a to 0.3 and testing gravity through the growth index 7 
to 3%. No external priors on the reconstruction parame- 
ters are necessary as they are selfcalibrated by the survey, 
most at the subpercent level. This also corresponds to 
subpercent calibration of the redshift space power spec- 
trum. 

We have tested the robustness of the conclusions by 
adding redshift evolution, which has little effect, vary- 
ing the number of wavenumber bins, and exploring the 
leverage from increasing the maximum wavenumber used, 
&max- Cosmological leverage plateaus by fc max = 0.4 — 
0.5ft./Mpc so the KLL form need only apply up to this 
scale. We made no assumptions about the departure 
from linearity, allowing the nonlinearity amplitude to 
float in a model independent manner in bins of k, but 
also then analyzed the impact of adopting a nonlinear 
prescription such as Halofit (or its revision). This im- 
proved the cosmology estimation and offers a promising 
sign to motivate continued development of cosmic emu- 
lators for the nonlinear power spectrum. 

Several areas exist for further development. The KLL 
form has been tested for dark matter, but [15] indicates 
it is successful for halos as well. Eventually this must be 
extended to galaxies, a major undertaking. On the posi- 
tive side, we achieved excellent results using reconstruc- 
tion starting from simple linear theory (Kaiser approx- 
imation); higher order perturbation theory approaches 
extend the range where reconstruction is milder. Uni- 
versality, i.e. cosmology dependence, of the reconstruc- 
tion is a major topic for future investigation, requiring 
large suites of cosmological simulations, again suited for 
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emulators. We have taken a first step toward address- 
ing this effect by exploring the influence of using Halofit 
and new Halofit cosmological dependences for the non- 
linearity. Simulations may also enable us to compress the 
information in bins down to a smaller set of parameters. 

Redshift space distortions provide a powerful tool for 
measuring the growth rate of cosmic structure, and de- 
livering insights on the competition between the gravita- 
tional laws driving clustering and accelerated expansion 
suppressing it. The results here give encouraging indica- 
tions, and quantitative measures, that theoretical anal- 
ysis can take into account robustly the nonlinear and 
velocity effects to extract accurate cosmology from the 
forthcoming large volume redshift surveys. 
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